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^ln  order  to  estimate  the  rate  of  vertical  transport  due  to  turbulence  In  the 
stratosphere.  It  is  necessary  to  know  the  degree  of  mixing  that  takes  place  In 
"Kelvln-Helmholt*  billow  events."  This  is  estimated  by  means  of  a  discrete 
model  and  the  results  are  compared  with  published  experimental  observations. 
From  these  considerations  It  is  concluded  that  a  very  large  degree  of  mixing 
probably  takes  place  In  such  events,  and  that,  therefore,  one  must  use  the  rela 
tion  between  the  bulk  vertical  eddy  diffuslvity  and  layer  diffuslvity  originally 
proposed  by  Rosenborg  and  Dewan  instead  of  the  simple  one  usually  employed.  ^ 
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Mixing  in  Billow  Turbulence  and 
Stratospheric  Eddy  Diffusion 


l.  INTRODUCTION 

The  value  of  the  vertical  eddy  diffusion  coefficient,  K^,  which  describes  the 
transport  effects  of  smull  scale  turbulence,  is  at  present  controversial  for  the  case 

of  the  stratosphere.  Two  relationships  have  boon  proposed  that  relate  bulk  14  to 

*!<  ® 

the  fraction  of  fluid  that  is  turbulent,  P  ,  as  measured  along  a  given  representative 
length  of  the  stratosphere.  One  of  these,  which  wbb  first  reported  by  Rosenberg 
and  Dewan, 1  is 

Ke  “  '  (1) 

— jr 

where  L  is  of  the  average  square  ofL,  and  L  is  defined  as  the  vertical  distance 
upwards  from  u  given  altitude  (assumed  locuted  somewhere  within  a  turbulent  layer) 
to  the  top  of  the  layer.  As  shown  in  Rosenberg  and  Dewun*  and  Dewan,**  the  manner 
in  which  L  Ib  calculated  yields  a  value  of  L,  that  would  bo  somewhat,  less  than  one  half 
the  thickness  of  the  turbulent  layers  on  average.  Rosenberg  and  Dewan*  gives 


(Received  for  publication  0  April  1070) 

1.  Rosenberg,  N.W.,  and  Dewan,  IS.  M,  (1075)  Stratospheric  Turbulence  and 

Vertical  ISffective  Diffusion  Coefficients,  Ar(.:i(L-TR-7b-0519,  AD  AO  19708. 

2.  Dewan,  E.  M.  (1079)  Estimate  of  Vertical  Eddy  Diffusion  Due  to  Turbulent 

Layers  in  the  Stratosphere,  ATGL-TlMO-oO^. 
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the  general  background  and  derivation  of  Eq.  (1)  while  Dewan  discusses  muny  de¬ 
tails  and  subtle  points  omitted  in  earlier  papers.  The  term  Atf  is  the  duration  of 
the  time  frame  of  the  model  developed  in  the  works  cited  above  and  its  value  is  es¬ 
timated  to  be  approximately  1800  sec.  The  other  relation,  which  is  used  by  Woods 
and  Wlloy^  and  Lilly4  is  given  by 


K 


e 


klp' 


(2) 


where  K.  designates  the  average  eddy  dlffusivity  within  the  turbulent  layers. 

It  was  shown  in  Dewan  that  Eqs,  (1)  and  (2)  represent  two  extreme  cases  that 
result  from  the  "vertical  stack"  model  for  vertical  transport  in  stratified  fluids  with 
intermittent  layers,  it  was  shown  there  that  Eq.  <  1  >  is  valid  when 

isi 

2 

where  AtM  is  the  mixing  time  which  is  of  order  Atf  (as  demonstrated  in  Dewan)  , 

Eq,  (2)  is  valid  when 

K,  At  «4,  IT.  (4) 

m  n* 

In  Dewan,  2  experimental  evidence  is  cited  in  support  of  the  condition  Eq.  (3). 

The  purpose  of  this  report  is  to  examine  this  question  theoretically.  It  will  be  con¬ 
cluded  that,  for  the  case  of  the  stratosphere,  the  model  for  mixing  developed  here 
supports  condition  Eq.  (3). 


2.  SIMPLE  MODEL  KOI!  MIXING  IN  A  TURBULENT  LAYER 


The  turbulence  to  be  found  in  the  stratosphere  is  primarily  due  to  the  Kelvin- 
Molmholtz  (K-H)  shear  Instability  in  the  context  of  stable  stratification.  The  pur¬ 
pose  of  this  section  is  to  derive  a  model  for  the  mixing  due  to  this  instability  in 
order  to  estimate  the  value  of  At^  or  its  equivalent,  so  thut  It  can  be  compared 
to[(4iT)/n2)j. 


3.  Woods,  J.D. ,  and  Wiley,  ILL.  (1072)  Hlllow  turbulence  and  ocean  micro- 

structure,  Deep  Sea  Research  and  Oceanic  Abst.  10)87-121, 

4.  Lilly,  D.K,,  Waco,  D,  1C. ,  and  Adelfung,  S.  1.  <  1975)  Stratospheric  mixing 

estimated  from  high-altitude  turbulence  measurements  by  using  energy 
budget  techniques,  The  Natural  Stratosphere  of  1074.  CHAP  Monograph  1, 
Final  Report,  IK )T-TST-75-5'l,  ppT'G-BV  to  G-'OO. 
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11  Definition  of  Kl 

is  of  course  not  constant  in  time  in  the  case  of  billow  turbulence.  In  con¬ 
trast  to  boundary  layer  turbulence,  a  billow  cannot  be  described  by  "steady-Btate" 
conditions,  Instead  it  consists  of  a  turbulent  breakdown  followed  by  a  decay. 

The  term  Atj^  is,  therefore,  somewhat  vague,  and  for  purposes  of  calculation 

we  sharpen  Its  meaning  by  defining 
at 

KL  AtM  "  /  KL  (t)  dtt  ,  W 

where  K“  Is  the  initial  eddy  dlffusivity  subsequent  to  breakdown,  and  (t)ls  the 
time  dependent  value  of  Henceforth  Eq.  (5)  or  it's  "finite  difference  equivalent" 
will  bo  ueed  in  place  of  "K£  At^. " 

(t)  is  defined  here  as 

Kl  (t)  ■  ui  ,  (6) 

where  u  is  the  turbulent  fluctuation  velocity,  that  is,  the  root  mean  square  deviation 

fromthe  mean  velocity,  and  l  is  a  representative  length  scale  for  the  energy  contain- 

&  6 

ing  eddya  (compare  Tonnekea  and  Lumley  ,  pp  44-47,  and  Paequlll  ),  We  shall 
assume  that  both  i  and  u  can  depend  upon  time,  The  decay  of  u  is  due  to  viscous 
dissipation  from  small  scale  motion  at  the  small  end  of  the  wave  number  spectrum 
cascade, 

To  define  l,  we  consider  two  plausible  possibilities)  The  first  is  (see  Tennekes 
and  Lumley8,  pp.  48  and  82) 

I  ■  5  .  (7) 

where  S  la  the  vertical  ahear  of  the  horizontal  mean  velocity.  The  other  is 

*  •  ft  .  (8) 

where  N  1b  the  buoyancy  frequency.  Eq.  (8)  follows  from  energy  considerations  to 
be  described  below,  In  the  development  of  our  model  below  we  will  choose  Eq.  (7) 
for  t .  This  will  be  done  In  order  to  get  a  lower  limit  estimate  of  the  degree  of  mixing, 

^This  Integral  converges  when  due  account  is  taken  of  viscosity. 

S.  Tennekes,  H. ,  and  Lumley,  J.  L„  (1912)  A  First  Course  In  Turbulence.  MIT 
Press,  Cambridge,  Massachusetts. 

8,  Pasquill,  F.  (1962)  Atmospheric  Diffusion,  Van  Nostrand  Co. ,  Ltd. 
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Next  we  consider  the  mixing  effects  of  Kj  (t)  upon  Itself  through  erosion  of  the 
mean  profiles  of  temperature  and  velocity  within  the  mixing  layer  as  the  turbulence 
mixes  material  within  the  layers;  tills  mixing  will  decrease  the  slope  of  mean 
quantities  with  respect  to  tile  vertical  coordinate  Z,  This  would,  by  definition, 
decrease  the  values  of  S  und  N.  This  in  turn  would  (from  Bqs,  (7)  and  (ti)  cause  * 
to  become  larger  and  In  turn,  from  Kq,  (6),  tond  to  increase  K,  (t>.  The  maximum 
value  of  h  (vertically)  Is  the  layer  thickness;  therefore*  if  at  the  start  Ji  has  this 
value,  thon  erosion  effects  will  not  lengthen  it. 

The  value  of  u  decreases  monotanlcally  as  turbulence  decays,  This  will  cause 
K,  (t)  to  tond  to  decrease  In  time,  as  will  be  shown  below;  however,  if  l  increases 

Lj  00 

or  remains  constant,  then  f  K,  (t)  dt  can  assume  much  larger  vulues  than 
•Hr  o  L 

(4  L  )/n  (unless,  us  willbe  shown,  L  Is  quite  small).  It  will  be  shown  that,  in 
this  way,  stratospheric  mixing  duo  to  billow  events  Is  probably  very  large  within 
the  layers  and  that  therefore  ICq ,  (3)  holds. 


2.2  Initial  Klueluallnn  Velocity 

The  Initial  fluctuation  or  rms  velocity,  Uj,  is  defined  as  the  velocity  of  turbu¬ 
lent  motion  immediately  after  the  K-H  instability  breaks  down,  It  will  dopend  upon 
the  Klchurdson  number,  which  Is  defined  as 


(0) 


where  g  Is  gravitational  acceleration,  <A  the  average  potential  temperature  over  a 
scale  as  large  us  the  entire  layer,  and  0  is  the  mean  potential  temperature  us  u 
function  of  height,  Z.  AU  is  the  difference  in  mean  horizontal  velocity  across  the 
layer,  and  AZ  Is  the  layer  thickness.  This  quantity,  as  defined  by  liq,  (0),  Is 
sometimes  culled  the  "layer  Klchurdson  number." 

To  estimate  u.  we  employ  u  simple  energy  budget  approach  that  is  based  upon 

,  1  7  U 

the  work  of  Busingor  and  Ludlum.  Figure  1  schematically  depicts  the  situation 
The  layer  of  thickness,  AZ,  has  a  shear,  S,  across  It,  Since  S  =  AU/AZ,  All  =SAZ, 
We  now  estimate  availability  of  kinetic  energy  of  this  configuration.  For  conven¬ 
ience  wu  define  8  U  a  AU/2  und  8  Z  s  AZ/2,  l.et  T7  be  the  mean  horizontal  velocity 
of  the  fluid  In  the  center  of  the  layer  just  prior  to  breakdown,  We  now  calculate  the 
kinetic  energy  made  available  by  the  exchange  of  the  four  parcels  us  Indicated  in 


7. 


B. 


Busingor,  J.A.  (1000)  On  tho  energy  supply  of  clear  ulr  turbulence,  in  Clear 
Air  Turbulence  and  Its  Petectlon,  edited  by  V.  11.  Puo  and  A.  Goldberg, 
Plenum’ Press,  New  S'orlt,  pp.  TOO- 10!!. 


l.udlam,  F.  il.  (  1007)  •  huractorlstlcs  of  billow  clouds  and  their  relation  to 
clear  air  turbulence,  Quart.  ,J,  Hoy,  Soc,  03;410-435. 

- - -  m 
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figure  1  (two  purculs  being  located  at  the  center,  and  the  other  two  on  opposite 
boundaries  of  the  layer).  The  resulting  available  fluctuation  kinetic  energy  (per 
unit  mass),  KE,  would  thus  be 


10  =  I 


(IT  +  Oil)  +  (IT  -  flU) 


2  _  2 

1  -  air 


=  2flu 


au.aS 


figure  1,  Parcel  Exchange  Diagram  for  Pillow 
Energy  Budget 


To  calculate  the  potential  energy  ehungo  (per  unit  mass),  PE,  caused  by  the 
procuss  of  displacing  the  four  pu reels  in  Figure  1,  we  assume  that  the  stable 
temperature  profile  is  not  ultcrcd  significantly  by  this  exchange.  We  ealculutu  the 
work  per  unit  mass  needed  to  move  one  of  these  parcels  u  distance  5  Z,  uguinst  the 
huoyanoe  force.  This  is  given  by 
555 


where 


which  Is  the  square  of  the  buoyancy  frequency,  N.  Four  parcels  thus  require  four 
times  the  work  in  Eq.  (U)i  therefore, 


iJ  lu!  J-C i-J  'J..  -  V.;  AuoAtiliiL-iLilf-i1 i -L  tl  I  -1  ii.e. 
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These  results  allow  an  estimate  of  Uj,  Subtracting  PE  from  KE  we  arrive  at 

the  total  remaining  available  kinetic  energy  that  can  give  rise  to  turbulent  velocity 

2 

fluctuations.  The  latter  will  be  3  (u^)  / 2  if  we  assume  isotropy.  ThuB 

4  ■  lift!!  (S2  -  N2)  ,  (14) 

hence 

h  -  h, 

U1  -SAzJ-j-A.  08) 

where 

Rj  ■  N2/S2  . 

2.3  Eitunition  of  Diwlpstion  Rate  and  Velocity 
Decay  Rate 

Having  an  estimate  for  Uj  we  now  estimate  u  as  a  function  of  time.  This  is 
needed,  of  course.  In  order  to  estimate  as  a  function  of  time  from  Eq.  (6),  As 
shown  by  Eqs.  (7)  and  (8),  there  is  an  Interrelation  between  4  and  u.  To  simplify 
this  complication  we  shall  divide  the  time  into  finite  epochs  during  which  4  can  be 
considered  ac  constant.  During  one  of  these  epochs  we  shall  also  assume  that  S 
und  N  are  constant.  The  duration  of  a  particular  epoch  will  be  taken  as  the  charac¬ 
teristic  time  for  the  fluctuation  motion  given  (from  dimensional  analysis)  by 

At  =  4/u  .  (16) 


where  4  and  u  are  the  values  at  the  start  of  the  epoch  under  consideration. 

Next  we  employ  the  well  known  relation  between  the  dissipation  rate,  e,  and 
the  velocity  and  length  scales  given  by 


From  tills,  one  can  calculate  decay  from 


d(2u2)  u3 

— as  =  '  t  • 


(18) 


Thus,  in  order  to  calculate  the  decrease  In  u  in  the  time  interval  At,  we  integrate 
the  separable  Eq.  (18)  and  obtain 
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where  u.  and  u.  are  the  initial  and  final  fluctuation  velocities  of  the  epoch  At  (com- 

*  Q  • 

pare  Dewan)  .  Using  Eq.  (18)  for  At  (u^  =  u  in  that  equation),  we  have 

uf  =  |  Uj  .  (20) 

We  shall  use  this  relation  below  for  an  iterated  sequence  of  epochs. 


2.4  Estimation  of  the  Erosion  of  S  and  N  and  the  Time 
Dependence  of  £ 

We  now  examine  the  erosion  of  mean  temperature  and  mean  horizontal  velocity 
profiles  caused  by  turbulent  mixing,  For  generality  we  shall  employ  the  symbol 
4(Z)  to  stand  for  either  the  mean  temperature  profile  6  (Z)  or  the  mean  horizontal 
velocity  profile  U(Z).  Figure  2  shows  our  assumed  linear  initial  profile  *(Z)  within 
the  layer.  Denoting  the  layer  thickness  by  A  ,  ^  letting  Z  =  0  at  the  bottom  of  the 
layer  and  letting  ♦  be  the  value  of  at  the  top  of  the  layer  with  ♦  »  Oat  the  bottom 
for  convenience,  we  have 


4  s= 


(21) 


Both  6  (Z)  and  U(Z)  within  the  layer  obey  the  diffusion  equation  during  mixing 
(assuming  constant  within  the  layer).  Thus, 

,2 


K, 


8“  y 
ezz 


a  * 

:TTT ' 


(22) 
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While  this  is  clearly  reasonable  for  6  (Z),  one  may  wish  to  consult  Batchelor 
186-189,  for  the  case  of  U(Z>,  The  initial  conditions  for  Eq.  (22)  in  the  region 
o  <  Z  <  A  is  given  by  Eq.  (21).  The  boundary  condition  will  be  taken  as 


pp. 


8  * 


T2 


Z--  0,  A 


0  . 


(23) 


t 

Formerly  denoted  by  AZ. 

9.  Dewan,  E.  M.  (1976)  Theoretical  Explanation  of  Spectral  Slopes  in  Stratospheric 
Turbulence  Data  and  Implications  for  Vertical  Transport.  AFGL-Tft-7(f-0247, 
AD'AO??63'07; - - — 

10.  Batchelor,  G.K.  (1967)  An  Introduction  to  Fluid  Dynamics.  Cambridge  Univ. 
Press. 
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Figure  2,  Initial  Density, 
Potential  Temperature,  or 
Velocity  Profile 


* 


A 


This  condition  represents  the  assumption  that  the  boundaries  of  the  layer  are  to  be 
regarded  as  perfectly  insulated.  The  observations  of  Woods  and  Wiley3  for  the 
ocean  lend  some  support  to  this  assumption,  and  it  would  presumably  hold  to  an 
even  greater  degree  for  the  atmosphere  (see  Section  4,2  below). 

The  standard  methods  for  solving  this  problem,  using  Fourier  series,  lead  to 


♦  (Z,t)  =  *o 


-  |K.  (n  tl 
e 


(24) 


where  g  (n)  *  -1  for  odd  values  of  n  and  zero  for  even  values.  Due  to  the  rapid. decay 
of  the  higher  harmonies  and  rapid  convergence  of  the  series,  we  need  consider  only 
the  dominant,  first  term.  Tims 


t)  ' 


j  (TT  /  /,  )2 


t) 


(25) 


Figure  3  shows  ty(Z,0)  us  well  as  'ii  (Z,  t)  during  diffusion.  Figure  3  is  only  a 

schematic,  not  a  mathematical  graph.  Next  we  wish  to  use  Kq,  (25)  to  determine 
2 

S  or  N  as  u  function  of  time.  Since  the  latter  quantities  are  proportional  to 
B't/S Z,  we  estimate  a  representative  value  of  aty/OZ  by  joining  the  end  points  of 
<MZ,  t)  at  Z  s  O  and  Z  =  A  as  indicated  in  Figure  3,  Thus, 


9*(t)  „  »(a  .  t)  -  jdo,  t) 

VZ  A 


(20) 
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Figure  3.  Profile  of  Mean 
Quantities  as  a  Function 
of  Time 


From  Eq.  (21)i 


ap-(4)V* 


-  KjHlMrt 


The  constant  (8/n  )  *  0.  811)  however,  we  shall  approximate  it  by  unity  in  order  to 
obtain  the  correct  Blopo  for  t  ■  O  at  the  beginning  of  the  first  epoch,  and  for  t  »  Att 

8*1  *0  -IKjUWM’Atl 

inrl  t*At  ~  t  e  •  (Z8) 


2  2  -  (Kr  or/  A)2  At) 

bP  (At)  »  N‘  (o)  e  L 


S  (At)  -  S  (o)  e 


-  (Kj^(n  /  A  )*  At) 


We  could  now  estimate  the  decay  of  N  or  S  from  epoch  to  epoch  on  the  buetB 
of  Eq.  (29).  Before  doing  so  we  first  derive  Eq.  (8)  from  energy  considerations. 
From  this  point  on,  u  will  refer  to  the  isotropic  velocity,  but  for  the  present  pur¬ 
pose  we  will  consider  it  to  be  in  the  vertical  direction.  Let  W  represent  the  work 
done  by  a  vertically  moving  parcel  as  it  goes  from  its  initial  equilibrium  position 
to  a  point  located  vertically  a  distance  J  away.  Then 


.niiuLEuto."* -.*tf 'ru.  L.  ..u,  >)  j 


(30) 


W  *  /  p  N2  ZdZ  , 


where  p  is  the  density  of  the  fluid.  Equating  W  to  the  kinetic  energy  of  the  parcel* 

2  |k, 2  1 2  * 

■£u_  ,  ,  (31) 

Hence* 


1  *  ft.  (32) 

The  important  point  to  notice  about  this  relation  it  thut  as  u  decays*  N  could  do- 
creaae  at  such  a  rate  that  4  remains  constant  or  even  increases  In  time.  ThiH 
consideration  applies  to  4  defined  by  Eq.  (71  as  well. 

2.3  A  Simple  Model  for  Mixing 

It  is  now  possible  to  set  up  a  flow  chart  for  a  discrete  time  or  "finite  difference" 
model  for  mixing  that  employs  all  of  the  above  observations*  so  as  to  give  as  a 
function  of  time. 

l.et  4j,  Uj,  Nj.  Sj,  and  Kj  be  the  Initial  values  of  the  respective  variables, 

We  have  hero  suppressed  the  subscript  L  in  K^.  Krom  previous  considerations 
we  havoi 

ix  »  Uj/N[  (33) 

from  Eq.  (32),  where  could  be  replaced  by  S^j 


Trom  Eq,  ( 15 ) ; 

Kt  -  ij  ut  (35) 

from  Eq.  (6),  and 

Atj  ■  41/u1  (36) 

from  Eq.  (18). 
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To  calculate  the  effect  of  erosion  of  N  or  S  in  one  time  Btep,  Eq.  (35)  gives 

n  n 


-  exp(-{KnAtn 


(37) 


where  the  vulue  n  -  1  designates  the  initial  step  and  the  factor  ■j  in  the  exponent 
would  be  absent  in  the  case  of  sn+i/Sn.  Subsequent  epochs  would  be  indicated  by 
successive  integral  values  of  n.  From  Eqs.  (20),  (8),  (8),  and  ( 16)  we  obtain  the 
following  relations  between  contiguous  epochsi 


n+ 1 


Jt 


"  hn 


“n+l^n+l 


n+1 

Kn+1  "  un+l  ^n+1 
Atn+1  “  inTl^Vl  * 


(30) 


Equation  (37)  completes  this  set.  Before  setting  up  the  flow  chart  for  this  model, 
it  is  useful  to  put  Eq.  (38)  into  dimensionless  form.  We  use  Nj  and  a  as  the  repre» 
sentative  scales  and  make  the  following  definitions) 


un 

Nn 

K 

At' 


un/(NjA) 

Nn/Nl 

Va 

W 


(30) 


where,  in  the  last  equation  of  this  set,  we  replaced  i  by  Its  value  given  In  Eq.  (38) 
(compare  with  Eq.  (30)). 


8.  STABILITY  AND  MIXING 

Our  objective  is  to  use  the  above  model  to  ascertain  when  the  criteria,  Eqs. 

(3)  and  (4),  are  valid.  When  this  is  done,  the  degree  of  molecular  mixing  will  be 
considered. 

3.1  Stability 

We  now  consider  the  stability  of  the  model  given  by  Eqs.  (38)  and  (37).  In 
order  to  arrive  at  a  conservative  estimate  of  what  layer  thickness  and  Kj  gives  rise 

1? 


) I—  III.  L. 


iIi.'l-'.-itV.iLa  rj  .  lIAtl  i .  ~i  (♦.'t.-iti  ■■  V  - 


to  effectively  total  mixing,  we  reconsider  the  model  In  view  of  the  ambiguity  of 
using  either  S  or  N  in  the  definition  of  i,  and  redesign  it  to  contain  the  moat  con¬ 
servative  aspects  of  both  approaches. 

For  example,  Eq.  (29)  gives  us 


S,.8,«p[-k1(»)Sm]. 


(40) 


(41) 


In  Eq.  (41)  the  erosion  could  be  significantly  faster  than  in  Eq.  (40)  due  to  the  fact 
that  the  absolute  value  of  the  negative  exponent  Is  twice  as  large.  To  be  conserva¬ 
tive  we  therefore  choose  the  milder  erosion  rate  given  by  Eq.  (40)  (N  is  an  inverse 
time  scale).  On  the  other  hand,  consider  the  case  where  R^  ■  0,  25  (the  threshold 
for  turbulence  in  moat  cases),  and  compare  the  values  of  i .,  the  initial  scale 
length,  for  the  two  definitions  of  1 1 


i 


l 


*7 


(42) 


ST 


(43) 


Thus,  in  the  case  of  Eq.  (43)  the  initial  mixing  would  be  less  than  that  for  the  case 
of  Eq.  (42)  in  view  of  the  fact  that  Kj  is  the  product  of  Uj  and  in  Eq.  (35).  Equa¬ 
tion  (43)  is,  therefore,  the  more  conservative  for  our  purposes.  Finally,  we  must 
choose  between  two  definitions  of  Atni  Atn  could  be  given  by  either  Sn  *  of 
We  seek  the  most  conservative  value  that  by  definition,  would  give  the  small¬ 
est  vulue  for  (Kn  At^)  in  the  exponent  of  Eq.  (37).  Since  S  >  N  for  <  1,  and  since 
wo  shall  only  bo  Interested  in  the  range  0  <  1(^  <  1  (stratified  turbulence  is  limited 
to  this  range),  it  follows  that  the  most  appropriate  value  of  Atn  Is 


At  *  1/S. 

n  n 


With  those  changes  we  now  define 
S 


S_ 


and 


*7 


(44) 


(45) 
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Experiments  with  the  hand  ealculutor  program  showed  thut  if  initial  Hj  were 
small  enough  to  cause  i  to  grow  from  epoch  to  epoch,  then  inevitably 

(X)  K'  At'  l  would  grow  indefinitely.  In  the  program,  l  was  replaced  by  the  upper 
n=l  n  "/ 

bound  A  when  the  calculated  vuluo  of  l  was  equal  to  or  larger  than  A  (that  is,  when 
l'  2  (1)  then  t '  5  1).  It  was  found  experimentally  thut  this  instability,  which  cuuaed 
J  to  grow,  occurred  when  Rj  *>  0.825.  Of  course,  it  is  also  useful  to  derive  this 
result  from  analysis, 

To  obtuln  this  result  algebraically,  we  first  prove  the  statement  that  if  l  re¬ 
mains  a  constant,  then 
n 

Limit  £  K' t'  -*«  .  (47) 

n^ce  n«i  n  n 

Using  Eq.  (44)  for  At,  Eq.  (38)  for  and  K  ,  Eq.  (39)  for  the  non-dimensional 
values,  Eq.  (47)  cun  bo  written 

t  ^  a  ■ 

n  u'i'  n  /  u '  \  2  « 

Lim  !  JC  B  Uni  ZJ  ( ■qr  1  K  S  (^n)  »  (48) 

n**«  I  r.-l  ”n  J  n^*®  [  n»l\  n/ J  n»l 


and  for  l  s  constunt,  the  sum  bucomes  equal  to  a  constant  X 
n 


(£,*) 


and,  hence,  is 


unbounded  as  n  becomes  arbitrarily  large.  To  find  the  crucial  value  of  Rj  we 
therefore  examine  successive  values  of  4(n.  In  particular,  we  consider  the  ratio 

of^n+l^n*' 

Using  l'  »  u'/S'andEq.  (40)  lor  (where  N  Is  replaced  by  S  )  so  that 
■  n  n  n  n  n  n 

(nohdimensionally) 


<  "p<- 


wo  find 


(Jr)’- 


(Figure  4  also  defines  these  terms)  whore  we  used  un+J  ■  (3/4  ur)  to  get  Eq.  (50). 
Initial  values  are 


(52) 


Sj  «  (H^'1  . 

To  explore  critical  behavior  we  set  j)2  =  d  r  V )  where  y  is  a  small  number. 

Then,  from  Eq.  (50),  for  n  =  1 


<i  +  y>-(|)  expdi2^2) 


log o  m)  (1  +y>] 


Next  we  shall  show  that  if  y  is  positive  then  l®  ut®°  Easter  than  unity. 

Using  Eq.  (50)  with  n  -  2, 


but,  since  ^l'1  +  und  usln®  .  we  have 

(  ^■)-(!)2«xp{n  +  y>1°«°  [(*)  (l  +  y)]j  • 


If  y  .  0,  (ig/i  t)  ■  1,  If  y  la  positive  it  is  easy  to  prove  (by  simple  expansions 
that  ig  <ig  >  1.  Furthermore,  In  a  similar  manner  it  can  be  easily  shown  that  if 
1 1  It  )  >  1  it  follows  that  l*ntl/*nl  >*•  s  By  mathematical  induction,  then,  we 
have  instability  if  (ig/ij)  >  1.  It  can  also  be  shown  that  if  d2/t\)  *  lf  the  vb1uo 
of  t  continues  to  decrease  and  one  has  stability.  Thus,  setting  y  »  0  in  Eq,  (54) 
and  determining  the  value  of  R^  (critical)  from  this  result  (that  is,  substituting  into 
Eq.  (51)  we  find  that  R,  -  0.825*  QED.  Since  R,  *  0.  25  is  assumed  here  to 

'chit  4 

be  necessary  for  the  initiation  of  turbulence,  it  follows  that  there  will  always  be  an 
instability  of  l  and,  hence,  total  mixing  according  to  this  model. 

*Thorpe1^  found  vertical  velocity  fluctuations  of  order  (^  )  at  the  start  of  billow 
turbulence  in  his  laboratory.  This  is  smaller  than  our  estimate  by  a  factor  of 

(yfs  )_1,  If  we  used  -  <^f  >\  "lA'  lnstettd  of  Eq'  (15)*  then  RlCR1T“  0,475 
instead  of  0.  825  and  our  conclusion  concerning  mixing  remains  unchanged. 

11.  Thorpe,  S,  A.  (1673)  Turbulence  in  stably  stratified  fluids!  A  review  of  labora 
tory  experiments.  Boundary  Layer  Met.  5164-119. 
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J  K(t)  (It  *  J  u  (t)  A  dt  =  A  J  Uj  dt  +  A  J  U|j  dt  . 


(64) 


Thu  a 


/  K(t)dtrJAJj  lOfi  u(%)  4  1  |  . 


(65) 


where  R  1  is  the  initial  Reynolds  number  [(U^  A  >/»a].  Thus,  Kq.  (65)  is  always 
finite  when  the  final  manner  of  decay  is  taken  into  account  (see  Eq.  (60)]. 

The  next  question^  is,  what  value  of  A  will  cuuse 

°°  ~~2 

/  K  (t)  dt  <  <  i-4-  (66) 

O  IT 


which  is  the  criterion  for  small  mixing  |Eq,  (4)],  In  Eq.  (66)  we  shall  ubb  the 
equality  for  definiteness.  In  that  CQBe,  UBing  Eq,  (65)  we  obtain  (using  L  -  ) 

3  [l0*e  (~ni)+  ‘]  ■  7  <87> 

Using  N  «  0,0225s'1,  v  «=  1.84X  10*4  m2/a  (standard  atmosphere  values  at 
20  km),  1^  *  0,  25  (hence  S  ■  0.  045s*  1),  and  using  Uj  a  Uj  in  Eq.  (34)  we  obtain 

A  (=  16.7  cm*  .  (68) 


This  is  about  threo  orders  of  magnitude  smaller  than  the  typical  vnluoH  of  A  re¬ 
ported  by  Rosenborg  and  Dowan  where  values  of  L  should  be  multiplied  by  2  or 
more  in  order  to  be  to  equul  A),  and  Anderson13,  Cudet13’ 14,  and  Uarat.  13,  13 
Therefore,  if  all  mixing  were  essentially  irreversible,  then,  from  the  above  results, 
one  could  conclude  thut  Eq.  (1)  and  not  Eq.  (2)  is  the  correct  one  to  use  in  the 
stratosphere,  We  must  therefore  examine  the  postulate  that  a  significant  amount 
of  molecular  mixing  takes  place.  Unless  this  latter  condition  is  met,  the  parcels 
of  fluid  would  migrate  back  to  their  stuble  altitudes.  ^ 


t 

'Note  that  if  one  were  to  arbitrarily  set  t^*,  then,  before, 

Unbounded,  f  finite  makes  liq.  (65)  finite. 

All  P"^ 

This  would  be  22  cm  If  u.  -■  ~r~  ttti — ~  used  (that  is, 

for  initial  vertical  velocity).  '  V  '  T 


«0 

J  K  (t)  dt  would  be 
o 


using  Thorpe's  value 


(Because  of  the  number  of  references  mentioned  in  the  above  text,  they  will  not  be 
listed  here.  Refer  to  Reference  Pago  29,  for  References  12  through  17.) 
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3.2  Degree  of  Molecular  Mixing 

We  now  wish  to  establish  whether  or  not  molecular  mixing  is  sufficient  to 

ensure  that  irreversible  mixing  tukea  place.  To  make  this  estimate  we  make  use 
1 7 

of  a  suggestion  of  Koop  to  compare  the  time  scale  of  the  K-H,  that  is,  the  "turn 
over  time,  "  T  y,  given  by 


to  the  molecular  mixing  time.  T  T  is  roughly  the  time  available  for  the  turbulence 
to  perform  u  significant  amount  of  mixing  uction. 

The  molecular  mixing  time  (or  diffusion  time),  r  is  estimated  by 


tD“  vD  /iD)  ' 


where  D  is  the  molecular  diffusion  constant  and  7)  ^  is  the  Kolmogorov  microscale 
for  the  cascade  (soe  Tennekes  and  LumleyR).  Irreversible  mixing  will  occur  when 
T  D  <  <  T  T  ‘  "I'°  [)  Wo  employ  the  viscous  microscale  given  by 

,,  1/4 

V  *(t)  •  (71) 


where  u  hero  must  bo  ropluccd  by  D  in  order  to  obtain?)  jy  D  Is  obtained  from  the 


Prandtl  number  P  ,  and  v  by 


«d*(£)  V3'4- 


P  *  0,  73  for  air  (Tennekes  and  Lumley, 5  and  for  water  P  ~  10  (Mill10); 

"  -4  2  " 

v  «  1.64  X  10  m  /s  for  the  lower  stratosphere  (20  km  — U.  S.  Standard  Atmos- 

lu  2  19 

phere)  ,  and  v  a  10  m  /s  for  water  (Hill  ),  From  Eq.  (7 1)  we  obtain  7)  ^  for 
the  utmosphere,  7)  jjA,  and  water,  ?)  Dw  respectively! 

18.  U.S,  Standard  Atmosphere,  1976.  NASA,  USAF,  NOAA  -  S/T  76-  11562, 

19,  Hill,  M,  N.  (196  2)  The  Soa,  Interscience  Pub. 


(74) 


V  DA 


1.  B4  X  10'3 


(75) 


Making  use  of  Table  1  In  Thorpe11  we  tabulate  here,  in  our  Table  1.  the  values  of 
t  that  are  typical  for  atmosphere,  ocean,  and  laboratory.  For  large  mixing, 

T  D  <  <  7  T‘  1  fives  the  values  of  all  the  relevant  quantities,  and  it  can  be 

seen  that  in  all  three  cases  T  D/r  <  <  1,  This  condition  holds  most  strongly  for 
the  case  of  the  atmosphere.  Thus,  it  seems  that  irreversible  mixing  is  assured 
in  the  case  of  the  atmosphere.  Thus,  Eq.  (3)  is  indeed  the  valid  criterion, 


Table  1.  Estimation  of  Irreversible  Mixing 
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4.  EXPERIMENTAL  MEASUREMENTS  IN  THE  LABORATORY, 
OCEAN,  AND  ATMOSPHERE 


4.1  Laboratory  ami  Oceanic  Observations 
20 

In  Woods  thero  is  reported  a  typical  profile  of  temperature  in  the  thermo- 
clinc  of  the  upper  oceun.  The  Step-liko  structures  seen  there  have  been  attributed 
to  mixing  events  thut  are  K-H  billows  of  the  type  dlscusBcd  here.  Businger7  des¬ 
cribes  the  mixing  and  spreading  process  of  K-H  ovents  and,  again,  total  mixing  or 
nearly  total  mixing  is  assumed  to  occur  within  the  layer,  In  view  of  the  analogy 
between  the  turbulence  in  the  upper  ocean  und  stratosphere,  it  follows  that  the  step- 
like  structures  seen  in  the  former  case  argue  for  the  existence  of  similur  structures 


in  tho  potential  temperature  profile  of  the  stratosphere. 

In  the  case  of  laboratory  measurements  there  have  been  no  observations  of  zero 
gradient  layers  in  density,  Thorpe11  has  consistently  found  an  approximately  linear 
gradient  across  the  layer  subsequent  to  breakdown  and  mixing,  The  same  effect 

20.  Woods,  J.D.  ( 1 08 B)  Wave-induced  shear  instability  in  the  summer  thermoclinei 
J,  Fluid,  Moch.  32:701-000. 
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occurs  in  the  observation  of  Koop.  Since,  in  fact,  there  seems  to  be  no  labora¬ 
tory  observations  of  the  "steplike"  structures  (which  in  turn  would  indicate  total 
mixing),  it  would  be  desirable  to  explain  the  discrepancy, 

Certain  differences  between  the  laboratory  conditions  and  those  found  In  the 
natural  environment  may  explain  the  apparent  contradiction.  For  example,  the 
experiments  of  Thorpe  and  Koop  involved  two  stably  superposed  fluids  that  were 
homogeneous  but  differed  In  density  across  a  very  narrow  interface.  The  homogen¬ 
eity  of  the  individual  fluids  would  allow  a  maximum  amount  of  layer  spreading  to 
take  place  after  the  onset  of  turbulence.  This  spreading  in  turn  would  cause  en¬ 
trainment,  and  the  latter  might  effectively  prevent  the  final  temperature  or  density 
profile  from  resembling  a  step.  The  step-like  quality  is  essential  in  our  model, 
since,  in  Section  2  it  was  assumed  that  the  boundaries  of  the  layer  acted  us  "insula¬ 
tors"  (see  Eq.  (23)].  It  is  therefore  unfortunate  that  a  mathematical  explanation 
for  the  discrepancy  between  laboratory  observations  and  the  oceanic  measurements 
is  not  yet  available. 

Table  2  based  on  Thorpe11  shows  the  values  of  three  other  parameters  relevant 
to  mixing,  namely,  R^,  the  Reynolds  number,  Ru,  the  Ruyleigh  number,  and  R^. 
TheBe  are  given  for  laboratory,  ocean,  and  atmosphere.  It  can  be  seen  that  Rg 
and  Ra  are  both  very  large  for  the  case  of  the  atmosphere  (in  comparison  to  the 
other  cases).  While  there  seems  to  be  no  clear-cut  theory  relating  these  parameters 
to  the  manner  that  mixing  lakes  place,  tile  re  Is  nevertheless  some  indication  that 
one  would  expect  thorough  mixing  in  the  atmosphere  (Thorpe11  related  Rg  to 
mixing). 


Table  2.  Parameters  Associated  With  Mixing 


4,2  Atmospheric  Observations 
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Perhaps  the  most  relevant  measurements  for  our  purposes  are  those  of  Mantis 
und  Pepping  These  consisted  of  radiosonde  observations  of  temperature  profiles 
in  the  stratosphere  and  troposphere.  Consider  the  following  quotations  from  their 
papers 


"Perhups  more  important. ....  is  the  observation  that  both  the 
stable  stratosphere  and  relatively  unstable  troposphere 
contain  numerous  shallow  layers  with  near  adiabatic  lapse 
rates, , . , ,  "  «• 

11 . characteristics  of  the  temperature  structure  lend  support 

to  the  hypothesis  that  there  is  a  change  in  the  turbulent  regime 
over  vertical  scales  of  order  100  meters.  The  deeper  unstable 
layers  have  these  dimensions  and  in  the  stratosphere  the 
unstable  layers  ure  frequently  capped  with  a  very  sharp 
Inversion. 11  ** 


These  remarks  ure  very  strongly  in  accord  with  the  model  proposed  here  and 
with  the  conclusion  that  large  amounts  of  mixing  take  place  in  K-H  events  in  the 
atmosphere.  The  step-like  structure  discussed  by  Woods  and  Wiley  has  also  been 
seen  in  the  stratosphere. 

Still  further  support  to  the  idea  of  large  mixing  in  an  "Insulated  layer"  will  be 

22  23 

found  in  the  radar  observations  of  Browning  and  Watkins  and  Atlas  et  al.  In 
both  cuses  a  double  layer  of  high  reflectivity  remained  for  a  long  time  after  the 
onset  of  turbulence.  This  can  be  considered  as  evidence  for  a  strong  inversion  at 
the  top  and  bottom  of  a  layer  that  has  undergone  complete  "insulated"  mixing  in  its 


I,  v;  interior.  The  layers  would  then  represent  the  edges  of  a  region  of  nearly  udiabatlc 

■  lapse  rate. 
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3.  CONCLUSION 

A  model  has  been  proposed  for  estimating  the  degree  of  mixing  due  to  K-U 
turbulence  events  in  the  stratosphere.  It  was  shown  that  this  model  indicates  that 


Hecall  that  complete  mixing  would  yield  adiabatic  lapse  rates. 

Complete  mixing  would  give  stoop  gradients  or  inversions  on  both  boundaries  of 
the  mixing  layer, 

21.  Mantis,  It.  T. ,  and  Pepin,  T.  J.  ( 197  l)  Vertical  temperature  structure  of  the 

free  atmosphere  at  mesoscale,  J,  Ueophs,  lies.  20i8621-8628. 
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22.  Browning,  K.A. ,  and  Watkins,  C,  D,  ( 197  0)  Observations  of  clear  air  turbulence 

by  high  power  radar,  Nature  227!260-263, 

23.  Atlas,  D.  ,  Metcalf,  ,1.1.,  lUehter?  J.  H. ,  and  dossard,  E.E,  (197  0)  The  birth 

of  'CAT1  and  microscale  turbulence,  J,  Atmos,  Scl.  27:903-013, 
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Urge  amounts  of  mixing  will  tatce  place  within  the  layer  and  that  it  la  irreversible 

due  to  molecular  diffusion.  This  conclusion  was  compared  to  laboratory  and 

environmental  measurements,  and  the  latter  were  found  to  give  further  support  to 

the  theoretical  prediction  of  large  mixing. 

In  conclusion^  Eq.  (2)  for  K  must  now  be  considered  invalid  and  in  its  place 

Eq.  (1)  should  be  used  to  estimate  K  in  the  stratosphere  (and  probably  also  in  the 

© 

ocean).  K„.  is,  of  course,  the  overall  vertical  eddy  diffusion  coefficient  for  small 
©  2 

scale  turbulence.  This  conclusion  casts  doubt  upon  the  previous  estimate  of  Lilly 

i  n 

at  al*  of  Km  m  0. 01  m*/s,  because  he  used  Eq.  (2)  to  perform  his  calculations. 

Note  added  in  proof 
Reference)  Bottom  p,  21 

The  value  Rj  ■  0. 625  was  obtained  by  using  the  initial  values  of  un  and  of 
eaoh  nth  epoch.  It  could  be  objected  that'  the  resulting  Kn  for  each  epoch  Is  over¬ 
estimated  and  that  therefore  the  calculation  overestimates  mixing.  It  is  thus 
Important  to  estimate  a  lower  bound  for  KR.  This  can  be  done  by  replacing  uR  and 
ln  by  the  lower  bounds  un  and  |  i for  edoh  epoch,  When  this  is  carried  out,  the 
resulting  critical  Rj  for  total  mixing  is  then  0.  689.  Our  conclusion  therefore 
remains  unaltered. 
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